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Context. Magnetohydrodynamic (MHD) waves are often reported in the solar atmosphere and usually classified as slow, fast, or 
Alfven. The possibility that these waves have mixed properties is often ignored. 
Aims. The goal of this work is to study and determine the nature of MHD kink waves. 

Methods. This is done by calculating the frequency, the damping rate and the eigenfunctions of MHD kink waves for three widely 
different MHD waves cases: a compressible pressure-less plasma, an incompressible plasma and a compressible plasma with non-zero 
plasma pressure which allows for MHD radiation. 

Results. In all three cases the frequency and the damping rate are for practical purposes the same as they differ at most by terms 
proportional to (k : R) 2 . In the magnetic flux tube the kink waves are in all three cases, to a high degree of accuracy incompressible 
waves with negligible pressure perturbations and with mainly horizontal motions. The main restoring force of kink waves in the 
magnetised flux tube is the magnetic tension force. The total pressure gradient force cannot be neglected except when the frequency 
of the kink wave is equal or slightly differs from the local Alfven frequency, i.e. in the resonant layer. 

Conclusions. Kink waves are very robust and do not care about the details of the MHD wave environment. The adjective fast is not 
the correct adjective to characterise kink waves. If an adjective is to be used it should be Alfvenic. However, it is better to realize that 
kink waves have mixed properties and cannot be put in one single box. 
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1. Introduction 

The last decade has seen an avalanche of observations of mag- 
netohydrodynamic (MHD) waves in the solar atmosphere. It is 
clear now that MHD waves are ubiquitous in the solar atmo- 
sphere. This has triggered new theoretical research for explain- 
ing and interpreting the observed properties. A special point of 
attention is whether these MHD waves are slow, fast or Alfven 
waves. Apparently, a large fraction of the solar MHD waves 
community favours very clear cut divisions and does not seem to 
appreciate the possibility of MHD waves with mixed properties. 
The trans verse oscillations observ ed in coronal loops (see for 
example lAschwanden et all Il999). often triggered by a nearby 
solar flare, are interpreted as fast kink MHD waves. A striking 
property of these transverse waves is their fast damping with 
damping times of the order of 3 - 5 periods. Resonant absorp- 
tion is up to today the only damping mechanism that offers a 
consistent explanation of this rapid damping. Resonant absorp- 
tion relies on the transfer of energy from a global MHD wave 
to local resonant Alfven waves. If this mechanism is indeed op- 
erational then this means that the observed transverse oscilla- 
tions have Alfvenic properties in at least part of the oscillating 
loop. The debate on the nature of MHD waves in the solar at- 
mosphere has gained new momentum when several groups, e.g. 



M. 



Goossens, 



e-mail: 



Send offprint requests to: 
Marcel . Goossens@wis . kuleuven . be 

* Postdoctoral Fellow of the National Fund for Scientific Research- 
Flanders (Belgium) (F.W.O.-Vlaanderen). 



|Pe Pontieu et all (|2007|) . lOkamoto et ail d2007l) . iTomczvk et all 
(|2007|) reported the detecti on of Alfven waves in HINO DE 
observations. Subsequently IVan Doorsselaere et al.l d2008l) ar- 
gued strongly against the possible presence of Alfven waves 
in the solar corona emphasising tha t Alfv en waves cannot be 
but torsional. IVan Doorsselaere et al.l d2008l) compared fast kink 
MHD waves to torsional Alfven waves and concluded that the 
HINODE observations have nothing to do with Alfven waves 
but can be explained in terms of fast kink MHD waves. 

This paper will not try to explain the HINODE observations. 
Its aim is to determine the nature of kink MHD waves on mag- 
netic flux tubes. We have no doubt about the explanation of 
transverse oscillation of coronal loops in terms of kink MHD 
waves. MHD waves with their azimuthal wave number equal 
to 1, i.e. m = 1, are the only motions that displace the axis of 
the loop and the loop as a whole. It is not clear to us on what 
arguments the use of the adjective fast is based. As far as we 
know there has not been any study of the forces that drive the 
kink waves in coronal loops. If these waves are fast, then the 
pressure gradient force should be, in general, the dominant force 
compared to the magnetic tension force. We have to admit that 
we also have used the adjective fast without a solid argument in 
favour of this classification. Our aim is to understand the spatial 
structure of the motions in the kink waves. 

An MHD wave on an axi-symmetric 1-D cylindrical plasma 
equilibrium is characterised by two wave numbers, the azimuthal 
wave number m, and the axial wave number k z . In addition 
modes can have a different number of nodes in the radial direc- 
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tion and this number can be used to further classify the modes. 
Hence, an MHD eigenmode is characterised by three numbers. 
The azimuthal wave number is an integer. The modes with m = 
are usually called sausage (slow and fast) or torsional (Alfven). 
The modes with m — 1 are named kink and the modes with 
m > 2 are flute modes. The axial wave number k z can be discre- 
tised as k : - re| with L the length of the loop and n = 1,2,... 
Depending on the dimensions of the equilibrium model there can 
be more than one radial eigenmode for a given couple (m, k z ). In 
what follows we shall study linear MHD waves that are superim- 
posed on a flux tube in static equilibrium with a straight and con- 
stant axial magnetic field. This equilibrium model contains the 
essential physics of the problem and allows a relatively straight- 
forward mathematical analysis of the MHD waves. MHD waves 
have been investigated in previous studies. However, these stud- 
ies almost exclusively focused on the frequencies of the MHD 
waves and in addition they were in m ost cases restricted to rea l 
frequencies. For example, the paper by Edwin & Rob erts! (1 19831) . 
which is often referred to in the solar MHD wave community, is 
limited to real frequencies and does not give any information on 
the eigenfunctions beyond the fact that they can be expressed 
in term s of Bessel fu nctions. Compl e x freq u encies were c onsid- 
ered by ISpruitl J 19821) and by ICallvl d!985h . ICallvl (1 19851) right- 
fully pointed out that essenti al physics is lost by restricting the 
analysis to real frequencies. Spruit] d!982l) is an exception to the 
rule in that he discussed the eigenfunctions. In short, the eigen- 
functions of MHD waves are not well documented even for uni- 
form equilibrium states and definitely not for non-uniform equi- 
librium states. 

In order to illustrate the nature of MHD kink waves we shall 
study them for three widely different MHD waves cases. The 
first case deals with compressible MHD waves of a pressure- 
less plasma on a high density flux tube embedded in a low den- 
sity magnetic plasma exterior. The assumption of a pressureless 
plasma removes the slow waves from the analysis and the veloc- 
ity and displacement have no axial component. The density be- 
ing higher in the flux tube than in the exterior means that accord- 
ing to common wisdom the MHD waves are propagating (body 
waves) in the flux tube and evanescent (surface wave) in the ex- 
terior. The second case deals with MHD waves of an incom- 
pressible plasma. The assumption of incompressibility removes 
the fast waves. The waves are evanescent (surface waves) both 
in the interior and exterior of the flux tube and this behaviour is 
independent of the density being higher or lower in the interior. 
In the third case we consider a situation where the compressible 
MHD waves in a non-zero beta plasma are propagating in the 
exterior. In that case the MHD waves are damped by MHD radi- 
ation, but the kink MHD wave hardly feels this wave leakage. 



2. Propagating kink MHD waves on dense 
pressureless flux tubes 

2.1. Equations for MHD waves for a pressureless plasma 
cylinder 

The equations for linear MHD waves on a 1 -dimensional 
pressureless cylinder with a straight fi eld can be obtaine d 
from the more gener al equations by e . g. lAppert et al.l (fl974l); 
Saku rai et aD dl991al) : iGoossens et all (119921) . iGoossens et al.l 
( 1995) by putting the local speed of sound Vs and the azimuthal 
component of the equilibrium ma gnetic field B w equal to 0. The 
resulting equations are (see e.g. iGoossensl 12008 ) 
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£, is the Lagrangian displacement and F is the Eulerian pertur- 
bation of total pressure. The coefficient functions D and C2 in 
(HJ are 



22/2 2 \ 
u — pv A u> \lo — co A \ , 



n 2/2 2 2 
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u> is the frequency and m and k z are the azimuthal and axial wave 
numbers, respectively. They define the dependence of the per- 
turbed quantities on time t and on the ignorable spatial variables 

(if, z) as 



exp[i(mtf + k z z - cot)]. 
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a»A and va are the local Alfven frequency and the local Alfven 
velocity respectively. They are defined as 
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For completeness we have also listed the expression for the local 
cusp frequency wc and the local speed of sound vs. For a pres- 
sureless plasma these two quantities vanish and as a consequence 
the slow modes are removed from the system since ij z = 0. 

The nature of an MHD wave is determined by the competi- 
tion of the restoring forces which are the force due to the total 
(gas plus magnetic) pressure gradient and the magnetic tension 
force. In ideal MHD we can obtain the following expression for 
the Lorentz force, the magnetic tension force IT and the magnetic 
pressure force V±P' in linear MHD waves on a background with 
a constant magnetic field 

i(VxB')xB 1- = -i^L +P o? A ^ 1, 
/ m , 2 \ 

n = - P u)\ \ r + ^ i y ) = -pw A £, 



dP' 

-V ± F = -—% 
dr 



-PX 



= -p (u> 2 - co A j l r + ijy 1^) 
= -p{o?-a> 2 A )t. 



(5) 



V ± is the gradient operator in horizontal planes perpendicular to 
the constant vertical magnetic field. The last equation in ((5} is 
obtained by using the second and third equations of (Q~|) to write 
the components of V±F in terms of the components of f . Note 
also that in the third line of ((5]) the first equality is general but 
the second equality is not. A similar thing happens between the 
fourth and fifth line. The second equality in both cases relies on 
the fact that g z — which is the case for a pressureless plasma. 
The important result from (0 is that the ratio of any of the two 
relevant components (the radial or azimuthal) of pressure force 
to the corresponding component of magnetic tension force is 
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In a non-uniform plasma with v 2 and hence co 2 K dependent on 
position, A also depends on position, meaning that the nature of 
the MHD wave changes according to the properties of the plasma 
it travels through. In a uniform plasma A is a constant and the 
nature of the wave does not change as it always sees the same 
environment. 

It is standard practice to rewrite the two first order ordinary 
differential equations of (fl~|i as a second order ordinary differen- 
tial equation for P': 



(y = k e r) in the exterior region. 



P[{r) = aJ m (x), 

h 



p((o z - aj A Ai ) 
P' e (r) = /3K ln (y), 

&Ar) = P - ^ 2 . K'Jy). 
p{co 2 - co 2 ) 



(14) 



tlp,X =\%-n*h\rl>, (7) 



dr 



where Y{io ) is the abbreviation for 



nco 2 ) 



k\ A(u> 2 ). 



(8) 



The prime denotes a derivative with respect to the argument x or 
y and a and ji are constants. Continuity of total pressure and the 
radial component of the Lagrangian displacement leads to the 
dispersion relation: 



J'Jx )K m (yo) , 
F = 1, 

J m (x )K^(y ) 

with the quantity F given by 



(15) 



F depends on or but also on the equilibrium through v\ and for 
a non-uniform plasma it is a function of position. The sign of 
F, i.e. of A, determines the local radial behaviour of the MHD 
wave. For F > the wave behaves locally in the radial direction 
as a propagating wave, for F < the the wave behaves locally as 
an evanescent wave (this is at least valid at large distances from 
the tube). 



2.2. Pressureless flux tubes with uniform density 
For a uniform plasma we can rewrite (0 as 



F = 



kj Pe(CQ 2 ~ (o\ e ) 
h pi{u) 2 - (J 2 ) ' 



(16) 



being xq = hR and yo = k e R. The dispersion relation ( fT5l ) can 
be solved numeri c ally. T his was done for real frequencies by 
Edwin & Roberts (1983) and for com plex frequencies by e.g. 
Spruitl (119821) and lcallvl dl985l 120031) . However, it is instruc- 
tive and also accurate to consider the so-called thin tube (TT) 
approximation (k-R <c 1). The Bessel functions J m (x) and K„,{y) 
in ( TTBT l are replaced with their first order asymptotic expansions. 
The dispersion relation (fT~5T > is reduced to 



d 2 P' 1 dP' 



dr 2 r dr 



+ -— - — -T(o/) P'=0. 



(9) 



Here T((o 2 ) is constant. We now specialise on uniform loops with 
a higher density p\ than the surrounding uniform plasma which 
has a constant density p e < p\. At the loop boundary r - R the 
density p changes discontinuously from its internal value p; to 
its external value p e . Since the magnetic field is uniform with 
the same strength everywhere, it follows that 



' Ai 
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,22,2 
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Note that the quantities T(w 2 ) and A(u> 2 ) are constant both in 
the interior and the exterior of the flux tube but change discon- 
tinuously at the loop boundary r — R. As our focus is on MHD 
waves that are propagating in the interior of the loop and evanes- 
cent outside the loop, we select the frequency so that 



Fi(w 2 )>0, F e (w 2 )<0, 
which means that 



(11) 



< a> < to a 



This allows us to define radial wave numbers h and k e as 



Ti(w 2 ) 



(O — CO 



Ai 



-r e (o; 2 ) 



or — a) 



Ae 



(12) 



(13) 



"Ae 



Equation (0 can then be solved in terms of Bessel functions 
J m (.x) (x = hr) in the internal part of the flux tube and K m (y) 



1 +F — =0. 

h 

The solution for the frequency is 

2 A^A, + P^je 2 
to = (O v , 

Pi + Pe 

and for the radial wave numbers k{ and k e 
kf = kt = k 2p ^^. 

~ Pi + Pe 



(17) 



(18) 



(19) 



(10) 

The right hand side of ( TT8l is almost invariably called the square 
of the kink frequency and denoted as a> 2 . In the thin tube approx- 
imation the frequ ency is indepen d ent of the wave number m > 1 
as already noted Gooss ens et al.l (Il992[) . Hence all flute modes 
with m > 2 have the same frequency as the kink mode with 
m = 1 . The radial wave numbers k\ and k e depend in a simple 
way on the density contrast. As the density contrast decreases, 
radial propagation is increasingly impeded and the nature of the 
MHD wave becomes gradually more Alfvenic, for p; = p e it is 
purely Alfvenic as the frequency of the global wave is equal to 
the local Alfven frequency everywhere. This is confirmed by the 
value of A 



Ai(w 2 ) = -A e (<o 2 ) = d—El, 

Pi +Pe 



(20) 



A(cj ) is constant in both the interior and exterior. It is positive 
in the interior meaning that the magnetic pressure force and the 
magnetic tension force act in the same direction. In the exterior 
its value is the exact opposite of that in the interior. The mag- 
netic pressure force and the magnetic tension force now oppose 
each other. Equation (f20b shows that the MHD waves are always 
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dominated by magnetic tension forces and that they are predom- 
inantly Alfvenic in nature. Take as an example a density contrast 
Pi/Pe = 3 then A; = 1 /2, A e = -1/2 so that the magnetic tension 
force is always twice as important as the pressure force. 
The TT approximations to the eigenfunctions are 



fr,i(r) 

R 
&,i(r) 

R 
P'(r) 



C, 
iC, 

C (k z Rj 



2 Pi 



V-6 = -C(LR) 



El L 

Pi + Pe R' 

2 P 



Pe r_ 

Pi + Pe R' 



(21) 



Note that when deriving (|2T1 > we have omitted terms of order 
(k z R) 2 and higher unless the terms of order (k z R) 2 are the first 
non-vanishing contribution to the expression under study. For 
example the expressions for £/ r ^{r)/R and %^i(r)/R mean that 
these two quantities are equal up to differences of order (k z R) 2 . 
The eigenfunctions are determined up to a multiplicative con- 
stant C which can be used to specify e.g. the radial displacement 
of the boundary of the loop. The radial and azimuthal compo- 
nents are tt/2 out of phase but they have equal magnitudes and 
in the loop they are constant. The wave is in the propagating 
domain in the internal part of the loop (k 2 > 0) but there are 
not any spatial variations. For realistic values of k z R = (ttR)/L 
the pressure perturbation and the divergence of the displacement 
field are zero for all practical purposes. The kink mode is to a 
high degree of accuracy an incompressible wave with very small 
magnetic pressure perturbations. Apart from f v the wave quanti- 
ties are continuous air - R. ^ varies in a discontinuous manner 
at r = R with opposite values at R K and R > . This discontinu- 
ous behaviour is due to the change of sign of the factor u> 2 - co 2 ^ 
when we move from the interior to the exterior of the loop. The 
behaviour of ^ creates strong shear layers which might undergo 
Kelvin -Helmholtz type instabilities as shown by iTerradas et alj 
(120081) . 

So far we have described the properties of modes that in- 
volve non-zero total pressure (P' + 0). However, if the medium 
is uniform the system of equations given by ([]]) also allows pure 
incompressible Alfven waves. They are only driven by magnetic 
tension, their eigenfrequency is simply <d = u>a and the total 
pressure, P', is equal to zero. To have such modes the displace- 
ment has to satisfy V-£ = 0. In cylindrical coordinates this means 
that, 



dr 



+ im^ v = 0. 



(22) 



If we prescribe the radial dependence of one of the components 
of the displacement, we can easily calculate the other component 
from the previous equation. The case m — is a particular solu- 
tion that represents torsional Alfven waves - 0, g v arbitrary), 
but equation ( l22l i can be solved for any azimuthal wavenum- 
ber i7i. Let us concentrate on m = 1 and the homogeneous loop 
model. Now we can have two different incompressible Alfven 
waves. One inside oscillating at the frequency o>Ai and another 
outside oscillating at u Ae . However, it is important to note that 
in such a configuration the radial displacement of the internal 
and external Alfven waves has to vanish at r = R (otherwise 
the continuity of the radial component is not guaranteed), i.e. 
the modes are localised in regions of constant Alfven frequency. 



This means that an internal incompressible Alfven wave is un- 
able to laterally displace the full tube (it cannot displace the tube 
boundary), although it is able to produce an incompressible mo- 
tion of the loop axis at its surroundings (for m = 1). Since these 
pure incompressible Alfven waves do not move the whole tube 
hereafter we will focus again on the kink solutions with P' + 0. 
Contrary to the incompressible Alfven waves these waves (with 
P' + 0) are able to connect the internal and external medium 
and to produce a coherent motion of the system because of their 
mixed nature. The fact that P' is small but different from zero 
plays a fundamental role in the mixed properties of these kink 
waves. 



2.3. Beyond the TT approximation for pressureless uniform 
flux tubes 

The analytic expressions ( |2TI ) (P' + 0) have been obtained in the 
limit k z R « 1. It is straightforward to solve the dispersion rela- 
tion (15[ and calculate the spatial solutions (TPfl i. This allows us 
to determine how the analytical expressions are modified by ef- 
fects due to a finite radius. In FigureQ]the eigenfunctions of three 
loops with different radii are represented. It is clear that the spa- 
tial profile is well described by the approximated solutions in the 
TT limit given by equations (fJTJ. The radial and azimuthal com- 
ponents are constant inside the loop, the azimuthal component 
has the expected jump at r — R, while the total pressure grows 
linearly with the radius. Increasing R results in an increase of the 
total pressure, and thus compressibility, since this magnitude is 
proportional to (k z R) 2 . Interestingly, for fat loops (see the case 
R/L — 0.1) the TT approximations of the eigenfunctions are still 
quite valid. An analysis of the forces (not shown here) indicates 
that, even for thick loops, the tension dominates over the mag- 
netic pressure gradient. 



2.4. Pressureless flux tubes with non-uniform density 

In this subsection we remove the discontinuous variation of den- 
sity from its internal value p\ to p e by a continuous variation 
in a non-uniform layer [R - 1/2, R + 1/2]. A fully non-uniform 
equilibrium state corresponds to / = 2R. When the jump in oja 
is replaced by a continuous variation of wa new physics is in- 
troduced in the system. The continuous variation of co A has the 
important effect that the kink MHD wave, which has its fre- 
quency in the Alfven continuum, interacts with local Alfven con- 
tinuum waves and gets damped. This resonant damping is trans- 
lated in a complex frequency (and complex eigenfunction). In 
the present paper damped global eigenmodes that are coupled 
to resonant Alfven waves in a non-uniform equilibrium state 
shall be computed by two methods. The first method is to use 
a numerical code that integrates the resistive MHD equations in 
the whole volume of the equilibrium state to determine a se- 
lected mode or part of the resistive spectrum of the system (see 



lected mode or part ot the resistive spectru m ot the system (see 
for example Van Doorsselaere et all [20041; lArregui et all [2005; 
Terradas et al The second method was introduced by 

iTirrv & Goossens! (119961) . It circumvents the numerical integra- 
tion of the non-ideal MHD equations and only requires numeri- 
cal integration (or closed analytical solutions) of the linear ideal 
MHD equations. The method relies on the fact that dissipation 
is important only in a narrow layer around the resonant point 
where the real part of the frequency of quasi-mode equals the 
local Alfven frequency. This makes it possible to obtain ana- 
lytical solutions to simplified versions of the linear dissipative 
MHD equations which accurately describe the linear motions 
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Fig. 1. Radial dependence of the normalised eigenfunctions, % r (real), 
g v (imaginary), and P' (real) calculated using d!4t after solving the dis- 
persion relation M5\ . The different lines represent three different loop 
radii, R/L = 0.01, 0.02, 0.1. The corresponding dimensionless eigenval- 
ues are u>L/v Ai = 3.847, 3.844, 3.803. The circles describe the incom- 
pressible eigenfunctions calculated using ( 140b after solving the disper- 
sion relation given by d4 1 b . The frequencies for the three loop radii are, 
in the incompressible limit, u>L/v A \ = 3.846, 3.842, 3.785. In these cal- 
culations pi/pe = 3 and k z = n/L. 



in the dissipative layer and in two overlap regions to the left 
and right of the dissipative layer. In these overlap regions both 
ideal MHD and dissipative MHD are valid. Asymptotic analy- 
sis of the analytical dissipative solutions allows to derive jump 
conditions that can be used to connect the solutions to the left 
and right of th e dissipative l a yer. Th e jump conditio n s wer e 
derived by e.g. ISakurai et ail (1 199 lab . IGoossens etail (1 19921) . 



iGoossens et alj d 1995b and IGoossens & Ruder marJ d 19951) for 
the driven problem and by iTirrv & Goossensl d 19961) for the 
eigenvalue problem. A schematic overview of the various 
regions involved i n this method is shown in Figure 1 of 
StenuitetalJC998). This method was used for computing eigen- 
modes of various no n-uniform plasma configu r ations by e.g. 
Tirrv et al ] (Il998allbh IStenuit et all (Il998l 1 19991) . lAndries et all 
(2000); Andries & Goossens (2001 ). A rela t ed m ethod was used 



bv lSakurai et alj (Il991bl) and lStenuit et al.l d 1995b for the com- 
putation of resonant Alfven waves in t he driven prob l em w ith a 
prescribed and real frequency and by Keppens et alj d 19941) for 
computing the multiple scattering and resonant absorption of p- 
modes by fibril sunspots. Comparison with results of fully dissi- 
pative computations show that the method is very accurate. 

A drastic variant of the method that avoids solving the non- 
ideal MHD equations uses the so-called thin boundary (TB) ap- 
proximation. In this lazy version the ideal MHD equations are 
not solved in the non-uniform plasma but the plasma is treated 
as if it were uniform all the way up to the dissipative layer. This 
is definitely a very strong assumption since the thickness of the 
dissipative layer is measured by the quantity 6a- 



6 A = 



1/3 



A A = 



dr 



(23) 



Here rj is the magnetic resistivity. If we denote / the typical 
length scale for the variations of the equilibrium quantities then 



OA 



(24) 



where R m is the magnetic Reynolds number (R m ~ rf l ). Strictly 
speaking the TB approximation assumes that the nonuniform 
layer and the dissipative layer coincide. It is common practice 
not to draw attention to this assumption but instead to refer to 
the TB approximation as the approximation t hat adopts l/R « 
1. Thi s TB approximation was first used by Hollweg & Yand 
( 1988 }). A discus sion of the TB approximation can be found in 
IGoossens! d2008h . 

In the TB approximation we need to add an additional term 
to the dispersion relation which takes into account the jump in 
the radial component across the resonant layer where the real 
part of the kink eigenmode is equal to the local Alfven frequency 
(jj = <jJa(^a)- fA is the resonant position which in the thin bound- 
ary approximation ta = R- The jump in £ r is 



[&] = -in 



m 2 /r 2 A 
P\Aa 



P', [P'] = 0. 



(25) 



The modified version of the ideal dispersion relation <JT5j is 

1. (26) 



J' m (x )K m (y ) . K m (y ) 
F 1 G 



J m (xoW m (y ) K' m (y ) 
F is given by (fTol l and G is defined as 

m 2 /r 2 A Pe (« 2 - <,) 



G — 71- 



P I A A 



(27) 



G contains the effect of the resonance. When we combine the 
thin tube (TT) approximation with the thin boundary (TB) ap- 
proximation, we can simplify the dispersion relation to 



k e k e R 
1 + F — - j G — = 0. 

ki m 



(28) 
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The zero order solution to d28l ). i.e. the solution when the effect 
of the resonance is not taken into account is of course ( TT~8T >. In 
order to take the effect of the resonance into account we write 



(29) 



and approximate a> 2 with oj^ + Ucj^j. The solution for the damp- 
ing decrement is 



n/2 m 



2 2 



«-<-.) 



2 \2 



Ae' 



cu k u>l R(p\+ p e ) 3 p(r A ) I A A (r A ) 



(30) 



Equati on (130b agrees with equation (77) of iGoossens et al.1 
dl992l) when that equation is corrected for a typo as the factor 
(ci) 2 Ki - o> 2 A ^ should be squared. This is surprising since that re- 
sult was obtained by IGoossens et al. I (119921) for surface waves 
in incompressible plasmas. In the same section of that paper it 
was noted that there is no distinction between compressible and 
incompressible plasmas for surface waves on thin tubes. 

Equation (f30b shows that the damping decrement depends 
linearly on m. Since we are mainly interested in m = 1 we shall 
specialise to that value in the remainder of this subsection. If the 
variation of cl> 2 a is solely due to the variation of density p as is the 
case here since we have considered a constant vertical magnetic 
field, equation ( 1301 can be rewritten as 



y_ _ nm (pi- p e ) 
aj k ~ 8 2? 



ft Ilk' 



(31) 



For a linear profile of density 



1 dr k 
so that 

_L 



Pi ~Pe 
I ' 



n I pi - p e 



8/?Pi+p e 



T 



4 1 Pi+p e 

n 2 l/R pi - p e 



(32) 



In ( l32l td is the damping time and T the period. Note that 
the result for ^ - of equation (l32l agrees with equation (79b) of 
IGoossens etlriTd 19921) . 

For a sinusoidal profile of density 



1 dr k 
so that 

y_ _ 

(L) k 



n Pi -p e 
2 / ' 



1 / p; - p e 



4 Pi + Pe 



td 

r 



2 1 Pi + p e 

nl/Rpi -p e ' 



(33) 



Here the results agre e with those obtained by 
iRuderman & Roberts! J2002). At this point we like to stress that 
the TTTB approximation turns out to be remarkably accurate 
far beyond its domain of applicability. This is cl early illustrated 
in a r ecent analytical seismological study by IGoossens et alj 
(2008) which complemented a fu lly numerical seismology 
investigation by Arregu ^et alj d2007l) . 

The eigenfunctions in the thin dissipative layer can 
be described by the functions F(r) and G(r) defined by 
IGoossens et al] (1 1995b fo r the driven problem an d the functions 



F(t) and G(t) defined by Rudem an et alj ( 19951) for the i ncom- 
pressible eigenvalue problem and by iTirrv & Goossens! (1996) 
for the compressible eigenvalue problem. In the dissipative layer 
the MHD kink waves are highly Alfve nic. This can be under - 
stood as follows. From the analysis by Sakur aTet alj d!991al) . 



IGoossens etail d 19951) and ITirrv & Goossens! ([1996) it follows 
that in the dissipative layer, the Eulerian perturbation of total 
pressure P' is constant and that | f r |«| t- v |. We do not have 
to worry about £. since it is zero for a pressureless plasma. In 
addition 



m P' F(t) 
r p 5 A | A A 



(34) 



Here t = (r - r A ) /6a is the stretched independent variable used 
in the dissipative layer which has a typical width [-55 A , 55a]- 
F(t) is a complex function and ^(r) does not suffer a discontin- 
uous jump as in the case of a uniform plasma. It is characterised 
by rapid spatial variation in the dissipative layer. From d34b it 
follows that the ratio of the <p component of magnetic tension to 
that of pressure gradient is 



on 



F(T) 



-p u\ %<p 

-i(m/r A )P' 



where we have used <jj 2 a / | A^ 
or bigger it follows that 



/ 



-i t F(j), 

OA 



~P W A 

-i(m/r A )P' 



(R m ) 1/3 » 1. 



(35) 

/. Since F(r) is of order unity 
(36) 



In the dissipative layer the magnetic tension force is far bigger 
than the pressure gradient force. Hence the MHD kink wave is 
highly Alfvenic in the dissipative layer. 

2.5. Beyond the TTTB approximation for pressureless 
non-uniform flux tubes 

In this subsection we go beyond the thin boundary approxima- 
tion and we consider thick non-uniform layers with l/R - 0.2 
and 0.4. The eigenfunctions are numerically calculated by solv- 
ing the full se t of lin ear, resistive MHD equatio ns described in 
lTerradasetaTld2006l) and using the PDE2D code (ISewelll[T995l) . 
Figure|2]displays the obtained results. It is clear that the two in- 
homogeneous solutions are almost identical to the homogeneous 
solution, except in the non-uniform layer, where large displace- 
ments are found. This is the location where the resonance takes 
place. 

In Figure [3] we see that both the radial and azimuthal com- 
ponents of the Lorentz force are dominated by magnetic ten- 
sion. The magnitude of the magnetic tension and pressure is of 
the same order (with the tension about twice as important as the 
magnetic pressure force such as corresponds to a density contrast 
of 3) except in the inhomogeneous layer where the tension in 
the azimuthal direction is clearly dominant, reflecting the strong 
Alfvenic nature of the solution in the dissipative layer. 



3. Incompressible MHD waves on flux tubes 

3. 1. Equations for incompressible MHD waves on a plasma 
cylinder 

The equations for in compressible MHD waves can be found in 
IGoossens et alJd!992l) . Incompressibility means that we take the 
limit vs — > °o and enforce V • £ — 0. Note also that in the incom- 
pressible case o>c = (jl>a- The relevant equations are 



P\o> 



dr 



m 



rP\ 
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Fig. 2. Radial dependence of the moduli of the eigenfunctions, 
and P' calculated numerically. The continuous lines represent the a- 
nalytical solutions for / = computed with d 14b (see also Figure QJ. 
Dashed and dotted lines represent the numerical computed eigenfunc- 
tions (R m = 10 6 ) for l/R = 0.2 and l/R = 0.4, respectively. The corres- 
ponding dimensionless eigenvalues are wrL/vai = 3.847, 3.799, 3.795 
and (DiL/v'Ai = 0, 0.092, 0. 192. In these calculations pjp s = 3, k z = n/L 
and R/L = 0.1. 



p(w 2 - aj 2 A )^ r 

p(a> 2 - w 2 a )%<p 
p(to 2 - u k )% z 



dP' 
dr 

—P', 
r 

ik 7 P'. 



(37) 



As before we concentrate on solutions with P' + (although 
now we assume incompressibility) and rewrite fl37l ) as a second 



Fig. 3. Top panel: Moduli of the radial forces. Middle panel: Moduli 
of the azimuthal forces. Bottom panel: Ratio of the modulus of the mag- 
netic tension to the magnetic pressure gradient in the radial and the azi- 
muthal direction. This case corresponds to the eigensolution displayed 
in Figure|2]with l/R = 0.4. The vertical dashed lines show the limits of 
the non-uniform layer. 



order ordinary differential equation for P' (eliminating the dis- 
placements): 




dP' 



(oj 2 - co 2 A ) 



dr 



(38) 



Here we do not need to worry about propagating and/or evanes- 
cent behaviour of the solutions and the local radial wave num- 
ber. The solutions are always evanescent or surface waves and 
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the local radial wave number is k,. Note that in our notation 

r(w 2 ) = -k\. 

3.2. Incompressible MHD waves on uniform flux tubes 
For a uniform plasma we can rewrite d38l l as 



d 2 P' 1 dP' { m 2 



dr 2 r dr 



(39) 



Equation ( |39l l can then be solved in terms of Bessel functions 
I m (x) (x = k z r) in the internal part of the flux tube and K m (x) in 
the exterior region. 



P[(r) = al m (x), 
&,i(r) = a 



k- 



Pi(to 2 - to 2 ) 



P' e (r) = [5K m {x\ 



k. 



p e (of - of) 



4W. 



K'Jx). 



(40) 



Continuity of total pressure and the radial component of the 
Lagrangian displacement leads to the dispersion relation: 



(41) 



(42) 



r m (x )K„,(x ) _ 

I m (x )K' m (x ) 
The incompressible version of F is 

_ p^of - af he ) 

Pi(W 2 - W 2 ;) ' 

Now xq = k z R. The dispersion relation ( |4TT > can be solved numer- 
ically. However, if we consider the TT approximation, the Bessel 
functions I,„(x) and K m (x) in < f4TT > are replaced with their first or- 
der asymptotic expansions, and the dispersion relation d4TT > is 
reduced to 

1 + F = 0. (43) 

The solution is again given by equation ( TT8l >. i.e. of = of^. 
The important point to note is that we get exactly the same 
expression for the frequency in the incompressible limit as in 
the pressureless limit. As far as the frequency is concerned the 
MHD waves with m > 1 do not make any distinction between 
a pressureless plasma and an incompressible plasma at least 
in the thin tube approximation. Here also A is given by d20i > 
Ai(a> 2 ) = -Ai(w 2 ) = S^-. This shows that the MHD waves 
are always dominated by magnetic tension forces and that they 
are predominantly Alfvenic in nature. Since the fast waves are 
absent as we have imposed the condition of incompressibility, 
the magnetic pressure force is associated with the slow wave be- 
haviour of the solution. The eigenfunctions are 



R 

R 

Zz,ir) 

R 
P[(r) 

(B 2 /p) 

&>,e(/?>) 



c, 

iC, 



iC(k z R)~, 



C (k 7 R) 



2 Pi ~Pe r_ 
Pi +p e /?' 



0, 

Hw(*<). 
-&(*<)■ 



As before we have omitted terms of order (k z R) 2 and higher un- 
less the terms of order (k : R) 2 are the first non-vanishing con- 
tribution to the expression under study. It is instructive to com- 
pare ( l44b with ( f2TT >. The expressions for % r j(r)/R, ^i(r)/R and 
Pj(r)/(B 2 /p) are exactly the same as in the pressureless case. 
There is now a small (of the order (k z R) axial component £. of 
the Lagrangian displacement. This axial component now makes 
V • £ exactly equal to zero while it is of order {k z R) 2 for a pres- 
sureless plasma. The azimuthal component ^(r) and the axial 
component % z j(r) are n/2 out of phase respect to £,-,i(r). £ r ,i(r) and 
£<pj(r) have equal magnitude and are constant in the flux tube. 
The axial component % z j(r) is of the order (k z R) and varies lin- 
early in the loop. In the thin tube approximation f-,;(r) is always 
small compared to £,,i(r) and ^,;(r). The dominant motion is in 
horizontal planes normal to the equilibrium magnetic field. The 
vertical displacement and velocity are small, of the order of k z at 
the tube boundary. As in the pressureless case ^ is discontinu- 
ous at r — R with opposite values at R K and /?>. This discontinu- 
ous behaviour is due to the change of sign of the factor cf - co 2 K . 
In the incompressible limit also f z is discontinuous at r — R with 
opposite values at R < and R>. This discontinuous behaviour is 
due to the change of sign of the factor of - of Q with cf^ 
the incompressible limit. 



cf A in 



3.3. Beyond the TT approximation for incompressible MHD 
waves on uniform flux tubes 

As in the compressible case the analytic expressions d44b have 
been obtained in the limit k z R « 1. It is easy to solve the dis- 
persion relation (PTTT i and determine the spatial solutions j40l . In 
Figure Q] the eigenfunctions of three loops with variable radii 
are represented with circles. Again it is clear that the spatial 
profile is well described by the approximated solutions in the 
TT limit given by equations d44b . As expected, the case with 
R/L = 0.1 shows the largest deviation from the TT approxi- 
mation. Figure Q] also allows us a direct comparison with the 
eigenfunctions in the compressible approximation. As the ratio 
R/L decreases the eigenfunctions of the compressible and in- 
compressible cases tend to be the same, in agreement with the 
analytical expressions given by (l44t and ( f2Tb . We arrive to the 
same conclusion for the forces, even for thick loops, the tension 
dominates over the magnetic pressure gradient. 

3.4. Incompressible MHD waves on non-uniform flux tubes 

We again remove the discontinuous variation of density from its 
internal value p; to p e by a continuous variation in a non-uniform. 
By doing so, we allow interaction of the global kink wave with 
local Alfven/slow continuum waves and the discontinuous be- 
haviour of % v and £, are replaced by singular behaviour in ideal 
MHD and by large but finite values in non-ideal MHD. In the 
thin boundary approximation we need to add an additional term 
to the dispersion relation which takes into account the jump in 
the radial component across the resonant layer where the real 
part of the kink eigenmode is equal to the local Alfven fre- 
quency. For the inco mpressible case the jump in £ r is (see e.g. 
iGoossens etaTlll992l) 



[&] = -in 
m = 0. 



m 2 lr\ + k\ 
pIAaI 



(45) 



(44) 



In the incompressible limit wa = wc and the Alfven resonance 
and slow resonance coincide. The jump in £ r due to the Alfven 
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resonance is given in d25l l. The jump in due to the slow reso- 
nance is 

,2 



[&] = -in 



P I Ac 



v l + v i 



[P'] = 0, A c 



d 
dr 



(a, 2 



3- 



(46) 



For the incompressible limit vi —> oo (146 1 becomes 



[&] = -in 



-,P\ 



IP'] 



P I A c | 
0, Ac = Aa 



d 
dr 



(a, 2 



(47) 



We see that equation ( |45T > is the combination of the jump due to 
the Alfven resonance (f25t and the jump due to the slow reso- 
nance (l47l i. The contribution to the jump in £ r due to the slow 
resonance is of order (k z R) 2 compared to that of the Alfven res- 
onance and can be neglected in our thin tube approximation. 
Note that (v|/(v? + v^)) 2 < 1 so that the slow resonance has its 
biggest effect for incompressible plasmas. Even in that case it is 
unimportant compared to the Alfven resonance. This is in agree- 
ment with a result obtained for MHD waves in prominences by 
ISoler et alJ d2009l) . For sentimental reasons we shall keep the 
contribution due to the slow resonance. 

The modified version of the ideal dispersion relation ( BTb is 



r m (xo)K m (x ) . K m (x ) 

t i G = 1. 

I m (x Q )K^(xo) K' m (x ) 

F is given by (fT~6b and G is defined as 

m 2/ri+fc 2 Pe(w 2 -w| e ) 



(48) 



(49) 



p | A A I k e 

When we combine the TT approximation with TB approxima- 
tion, the dispersion relation is reduced to 

k,R 

1+F-iG — = 0. (50) 
m 

The zero order solution to (l50l without taking into account the 
effect of the resonance is of course (fT8l . The effect of the reso- 
nance is contained in G. In order to take that effect into account 
we proceed as before 

a) — wr + iy, wr = (tit, (51) 

and approximate aj 2 with u^ + licj^y. The solution for the damp- 
ing decrement is 

y _ 7t/2 (m + (kzRf \ P 2 pl Ki-Q 2 
Wk oj\ \R mR J (pi + p e ) 3 p(r A ) | A A (r A ) |" 

If the variation of to 2 k is solely due to the variation of density p 
as is the case here since the vertical magnetic field is constant, 
the equation can be rewritten as 



y_ 



n\m (k z R) 2 \ (A-pe) 2 1 



R mR j p i+ p e | f \ rx 



(53) 



From here on we shall specialise to m = I. For a linear profile of 
density 

JL = JL( l + {kzR f)LPlzPi, 

co k 8 V - 'Rpi+pe' 

IE - 4 l A + p. 

T ~ 7T 2 l/Rp i -p l 



^(i-(M) 2 ). 



For a sinusoidal profile of density 

OJk 4 V " 'Rpi+pe 

T£ _ 2 1 pi + p e 

T 



■ (l - (k z R?) ■ 



(55) 



It I/ R pi - p e 

For all practical purposes we can neglect the contribution pro- 
portional to (k z R) 2 and conclude that the damping due to reso- 
nant absorption of the kink mode in an incompressible plasma is 
the same as that in a pressureless plasma (see equations fl32| l and 
(l33ll). If we forget about differences proportional to (k z R) 2 then 
the conclusion is that kink MHD waves in pressureless plasmas 
and incompressible plasmas are the same. In view of that conclu- 
sion it is difficult to understand why a kink mode can be called 
fast as fast waves are absent from incompressible plasmas. 

The eigenfunctions in the thin dissipative layer can be de- 
scribed b y the functions F(t) a nd G(t) which were first intro- 
duced by Rud eman et al.l d 1995b for non-stationary incompress- 
ible resonant Alfven waves in planar plasmas. The conclusion is 
the same as in the previous section. The kink MHD waves are 
highly Alfvenic in the dissipative layer. 

4. MHD kink waves in the presence of MHD 
radiation 

4.1. Equations for compressible MHD waves on a non-zero 
beta plasma cylinder 

So far we have seen that kink MHD waves in the thin tube 
approximation do not care about propagating (body wave) or 
evanescent (surface wave) behaviour in the internal part of the 
flux tube. The behaviour in the exterior plasma was until now 
evanescent. Here we take the next step and consider leakage of 
energy due to MHD radiation. MHD radiation causes the fre- 
quencies to be complex even in absence of resonant damping. 
MHD waves in the prese nce of MHD radiation were studied 
for unifor m flux t ubes bvlSpruitl dl982l) in the TT approxima- 
tion and bvlCallv| | 19851 f2003 ) for arb i trary v alues of the radius. 
IStenuit et al.1 (Il998h and IStenuit etall d 19991) determined MHD 
waves undergoing resonant absorption and/or leakage for pho- 
tospheric flu x tubes embedded in a non-magnetic surrounding. 
IStenuit et al.l (Q 999) pointed out which Hankel function to use 
for leaky and non-leaky waves. We use the equations for linear 
MHD waves on a 1 -dimensional cylinder with a straight field. 
Effects due to plasma pressure and compressibility are taken into 
account. The equations are 



D- 



p\o) 



dr 

,x dP' 



in i\ im , 

p(</-<ayf, = —P\ 



pyoj 



4) ^ 



ik- 



P'. 



-co 2 P' 



The coefficient functions D and C2 are now 

D = p(v 2 +v 2 A )(co 2 -co 2 A )(oj 2 -aj 2 ), 

(54) C 2 = a/-(v 2 + v A )(^-4)(^ + 4 



(56) 



(57) 
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As before we rewrite the two first order differential equations of 
d56l l as a second order ordinary differential equation for P'\ 



dp' 



[p(cj 2 -cj 2 a ) 



dr 



— -Y(u 2 )\rP\ (58) 



(59) 



where T(oj) is now defined as 

r(w ) = — — — . 

(v 2 s + vi)(^-4) 

We have solved the set of equations (|56l l under general con- 
ditions allo wing f o r non- z ero plasma p ressure and compress- 
ibilit y, see ISpruitl (fl982l) : ICallvl (1 19851) for uniform plasmas 
and lGoossens & Hollwedd 19931) for nonuniform plasmas. Here 
we present the results for a pressureless plasma with Vs = 0. 
Equations (|56] l are then reduced to equation ([TJ and equation 
d59l l is reduced to equation ((8). On one hand, sine we aim to 
study MHD waves that are propagating in the external medium 
we require that oj 2 - tJ A > 0. On the other hand, we want to 
show that kink waves do not need trapping and since we want 
to have resonant absorption present in the model we require that 
lo 2 - (jj 2 Ki < 0. In case of a constant magnetic field then the in- 
equalities imply that pi < p e so that the flux tube is underdense. 
Note that when we allow vs + and/or the equilibrium magnetic 
field to be non-constant then MHD radiation does not necessar- 
ily require an underdense loop. 



4.2. MHD waves on uniform flux tubes 



where F is defined by ( fTST l. Now xo = k,R and yo = k e R. The dis- 
persion rejation d62l can be solved numeric ally. This was done 
by e.g. ISpruitl (1 19821) and lCallvl (Il985ll2003l) . Using the TT ap- 
proximation the dispersion relation d62l is reduced to 



i + f|{i + ^} = o. 



(63) 



When we neglect in a zeroth order approximation the effect 
of MHD radiation the solution to equation (1631 is again equation 
(fT8l for the square of the frequency and equation ( fT9] > for the 
radial wave numbers. 

We can rewrite equation d63l correct up to second order in 
(keR) as 



k e n r. 
k[ 2 
The solution is 



(64) 



7 



PiPe 



— = 7^) 2 - -r 

OA 4 (Pi+Pe) 



(65) 



Equation ( I65l l then takes the simple form 

x = _jw£l^. (66) 

w k 8 (pi + Pe) 

The ratio of the force due to the pressure gradient to the magnetic 
tension force is here also given by equation ((6]). 

The TT approximations to the eigenfunctions are again given 
by equations ( |2T1 ). If we allow vs + then we find a non-zero £. 
with 



For a uniform plasma without gas pressure we recover equation 
©. Now Tiioj 2 ) < and F e (w 2 ) > and w 2 e < u 2 < gj*. The 
radial wave numbers are now defined as (see equation (1 131 )) 



R 



(k z R) -. 
R 



(67) 



2 2 2 2 

. - 0} — UJ.. - . to — OJ\ 

k 2 = -ri(co 2 ) = 5-A ^ e 2 = r e ( w 2 ) = — r-t*.. (60) 



"Ae 



The solutions to equation (O are now 



P[(r) = al m (x), 
&,i(r) = a 



1 



Kif) = /3H^(y) 



p(oj 2 - oj 2 Ai ) 



&, e (r) = P 



1 



P^ 2 ~ W Ae) 



h KM), 



kH'V(y). 



(61) 



The formulation with the Hankel functions is convenient as it 
enables us to distinguish between incoming and outgoing waves. 
With the classic convention that Re(oj) > 0, and the time depen- 
dency exp(-/«f), the function H$ corresponds to an outgoing 

(2) 

wave and H;„' to an incoming wave. Hence we drop a possible 
contribution due to H*^ since we are dealing with the eigenvalue 
problem and not with the driven problem where the incoming 
wave is prescribed. For the solutio ns interior in the tube we have 
taken the Bessel function I m as in lGoossens et al.l d!992l) . 

Continuity of total pressure and the radial component of the 
Lagrangian displacement leads to the dispersion relation: 



p r m (xo)H£>(y ) 
7m = ' 

/ m (*0)#m (Jo) 



(62) 



4.3. MHD radiating waves on non-uniform flux tubes 

Once more we remove the discontinuous variation of density. 
Now the jump in £ f is given by d25l ). and the modified version of 
the ideal dispersion relation d62l is 



(68) 



The TT approximation to this equation (for m = 1) is 



1 +F 



k e (. .TT 
— < 1 + l- 

h \ 2 



(k. : R) 2 \ - /G(*e/?){l + £(k e R) 2 



0. 



(69) 



In the third term of the left hand side of the previous equation 
we can drop i^(k e R) 2 for two reasons. First it produces a term of 
order (A; e R) 3 , second it leads to a (small) change in the real part 
of the frequency wr which we approximate by o)k. We can then 
rewrite ( f69b correct up to second order in (k e R) as 



(70) 



1 + F -T = M^efl) 2 + iG(k e R). 
k[ 2 

The solution of ( |70l is 

WR = OJ k , 

y n/2 1 pfpl (o^-cojj 2 



w k u>l R (Pi +p e ) 3 p(^a) I A A 



(Pi +Pe) 2 



(or 



(71) 
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In case of a constant magnetic field so that the variation of the 
local Alfven frequency is solely due to a variation of density (fTTT i 
can be further simplified to 



JL B JLLezz* - fW<fr-*> < o. 

&Rpi+p e 8 (pj + p e ) 2 



(72) 



Equation d72l is derived for a linear variation of density. For a 
sinusoidal variation the factor 7r/8 is to be replaced with 1/4. 
The absolute value of the ratio of the decrement due to resonant 
absorption (first term of equation ( l72b ) and that due to MHD 
radiation (second term of equation d72"l )) is 



l/R p e + pj 
(k z R) 2 p e - Pi 



(73) 



This clearly indicates that damping due to resonant absorption 
dominates over that due to MHD radiation (since k z R « 1). For 
example, for a tube with p e /pi = 3 and R/L = 0.01 then a very 
tiny layer of l/R = 5 x 10~ 4 is enough for resonant absorption to 
dominate over radiation. 

The remarkable result is that the frequency of the kink wave 
and its damping due to resonant absorption found for a com- 
pressible plasma with a non-zero plasma pressure differ by terms 
of order {k z R) 2 , even when we allow MHD radiation. If we ne- 
glect contributions proportional to (k z R) 2 then the simple con- 
clusion is that the frequency of the kink wave and its damping 
due to resonant absorption are the same in the three cases that 
we have considered. 

Again, as in the two previous cases, the eigenfunctions in the 
thin dissipative layer can be described by the functions F(t) and 
G(r). Again the conclusion is that kink MHD waves are highly 
Alfvenic in the dissipative layer. 

5. Conclusion 

This paper has examined the nature of MHD kink waves. This 
was done by determining the frequency, the damping rate and 
in particular the eigenfunctions of MHD kink waves for three 
widely different MHD waves cases: a compressible pressureless 
plasma, an incompressible plasma and a compressible plasma 
with non-zero plasma pressure which allows for MHD radiation. 
The overall conclusion is that kink waves are very robust and do 
not care about the details of the MHD wave environment. In all 
three cases the frequency and the damping rate are for practi- 
cal purposes the same as they differ at most by terms propor- 
tional to (k,R) 2 . In the magnetic flux tube the kink waves are 
in all three cases, to a high degree of accuracy incompressible 
waves with negligible pressure perturbations and with mainly 
horizontal motions. The main restoring force of kink waves in 
the magnetised flux tube is the magnetic tension force. The gra- 
dient pressure force cannot be neglected except when the fre- 
quency of the kink wave is equal or slightly differs from the local 
Alfven frequency, i.e. in the resonant layer. The adjective fast is 
not the correct adjective to characterise kink waves. If an adjec- 
tive is to be used it should be Alfvenic. However, it is better to 
realize that kink waves have mixed properties and cannot be put 
in one single box. 
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